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Abstract: We show that the single emitter linewidth underlying a broad- 
ened ensemble emission spectrum can be extracted from correlations among 
the stochastic intensity fluctuations in the ensemble spectrum. Spectral cor- 
relations can be observed at high temporal and spectral resolutions with a 
cross-correlated pair of avalanche photodiodes placed at the outputs of a 
scanning Michelson interferometer. As illustrated with simulations in con- 
junction with Fluorescence Correlation Spectroscopy, our approach over- 
comes ensemble and temporal inhomogeneous broadening to provide single 
emitter linewidths, even for emitters under weak, continuous, broadband ex- 
citation. 

OCIS codes: (030.5290) Photon statistics ; (300.6280) Spectroscopy, fluorescence and lumi- 
nescence. 
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1. Introduction 

Single emitters often display dynamic and complex behaviors entirely masked by ensemble 
measurements. Spectroscopy on these emitters has reached impressive heights. There are sys- 
tems, however, where it is not feasible or desirable to separate the individual emitter from the 
ensemble. When confronted with the task of isolating individual properties from large popula- 
tions, spectroscopy offers a variety of dedicated responses such as Doppler-free, hole-burning 
and photon-echo spectroscopies IT] to resolve ensemble and temporal averaging effects encoun- 
tered in inhomogeneously and homogeneously broadened samples. These diverse and powerful 
techniques share a common trait; they all rely on the nonlinear optical properties of the emitters. 
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Fig. 1. Our approach uses an interferometer to convert fast spectral fluctuations into in- 
tensity fluctuations. The cross-correlation of the fluctuations in the interferometer outputs 
reveals single emitter spectral dynamics with high temporal and spectral resolution. 



Hence, they tend to perform poorly on emitters with a small or vanishing optical nonlinearity 
and questionably on samples too delicate to handle the high excitation power required for non- 
linear optics. 

Surprisingly, methods to extract the linewidth of a single emitter from a broadened ensemble 
spectrum under the more gentle conditions of linear excitation remain more elusive. In prin- 
ciple, the single emitter linewidth can be determined from correlations among the stochastic 
fluctuations in the ensemble emission spectrum |2|. Progress in this direction was demonstrated 
in previous works investigating the autocorrelation of the broad spectrum of disordered nanos- 
tructures Due to the long duration necessary to record a high resolution spectrum, that 

approach cannot resolve the fast temporal broadening effects found in most inhomogeneous 
samples, severely limiting the ability to measure an underlying single emitter linewidth. 

In this Letter, we describe an experimental method revealing spectral correlations of a single 
emitter with high spectral and temporal resolution, despite the single emitter spectrum being 
obscured by an ensemble emission spectrum. The approach - shown in Figure 1 and reminis- 
cent of our previous work on Photon Correlation Fourier Spectroscopy |5] - consists of using 
a scanning Michelson interferometer to turn spectral correlations in the broadened spectrum 
into intensity correlations recorded by a Hanbury Brown Twiss detection setup. Previously, we 
applied this method to an isolated single emitter. Here, we expand and generalize to obtain the 
same dynamic single emitter spectral information from an ensemble of emitters. The method is 
investigated theoretically and demonstrated with numerical simulations illustrating its signifi- 
cance in conjunction with Fluorescence Correlation Spectroscopy (FCS). 

2. Intensity correlations in an inhomogeneous spectrum 

We introduce here the various quantities necessary to describe intensity correlations in an in- 
homogeneous spectrum and then describe a potential experimental setup dedicated to their 
measurement. 

2.1. Theoretical description 

We consider a collection of nearly identical emitters embedded in a dynamic and spa- 
tially inhomogeneous environment. Due to slight variations in the structure of the emitters 
(e.g. in their shape or composition) and spatial inhomogeneity (as caused, for example, by 
nanoscale disorder for emitters in a condensed medium), the spectra of the single emitters do 
not collapse on a single narrow, homogeneous stationary spectral line, but instead disperse 
their averaged frequencies coi , 0),, OJat over a spectral width A around some frequency 
(Oq (Figure 2a). Temporal inhomogeneities additionally cause the lineshapes of the emitters, 
s\{(0,t),...,Si{(0,t),...,SM[(0,t), to be explicitly dependent on time f, each emission fine un- 




Fig. 2. Single-emitter and ensemble spectral correlations in the fluorescence of a population 
of emitters, a) We first consider an ensemble of emitters with the same narrow line shape 
but different center frequencies. The ensemble spectrum would appear as a broad Gaus- 
sian (5(a), f)). b) With our setup, intensity correlations from the emission of single emitters 
flowing under a microscope objective will contribute to the single emitter spectral correla- 
tion, ps™g''=(^, t), if the two photons correlated are from the same emitter. Correlations of 
photons originating from different emitters will provide the ensemble spectral correlation. 



dergoing independent, identically distributed temporal stochastic spectral fluctuations over a 
mean-squared range = ([i,(a),f) — fi),]^) centered around the single emitter's average transi- 
tion frequency 0), (where (.) denotes the average over many independent observations). At any 
time t, each emitter contributes to the ensemble emission with an intensity /,(?) / Si{(0,t)dCO 
and a normalized lineshape .?,■(©, f) = Si{(0,t)/Ii{t) (i.e. /.f/(a),f)da) = 1). 

Stochastic fluctuations in the total intensity of the spectrum I{t) = L/Lo^'(0 ^^e tradition- 
ally approached through the second-order (intensity) correlation function ^'^'(t) ~ {I{t)I{t + 
t)) / {I{t))^. Similarly, fluctuations in the ensemble spectrum 5((B,f) = L/1o'^'(®jO be an- 
alyzed through the spectral correlation function t) : 

P(C,t) = ( j 5(a),f)5(co + C,f + T)dco). (1) 

Qualitatively, the spectral correlation function t) scales as the probability to measure a fre- 
quency difference between two photons separated by a time interval T. Important properties 
of f (i^, t) appear when splitting up EqlUinto two distinct components as : 

P(C,T) =P™-'(C,t)+P™'s'"(C,t) (2) 

where ^^'(Ct) = Zi^j j {si{(o,t)sj{(o + C,,t + and P™gi'=(C,T) = A^(/i;(a),f)i/(« + 

1^,? + T)d(B). Since emitters have statistically independent fluctuations, the component 
P™'(C,t) rewrites as f"="'(C,T) = j {si{co,t)){sj{(0 + C,,t))d(0, i.e. ^""'(Ct) is indepen- 
dent of T and reduces to the autocorrelation of the average ensemble spectrum. Similarly, 
psingie^^^ t) reduces to the autocorrelation of the time-averaged single emitter spectrum as long 
as the single emitter spectra Si{(i),t) and Si{co,t + t) are uncorrelated, i.e. for large values of 
T. But on shorter timescales T ^ 0, temporal inhomogeneous broadening is suppressed as no 
spectral fluctuations have time to occur, and P^'^^^'^{^,z = 0) coincides with the autocorrela- 
tion of the single emitter spectrum. The difference between P™'s''^(i^,T = 0) and is 
illustrated in Figure 2b. 

Information on the linewidth of the single emitter hence appears systematically encoded 
into the spectral correlation function P{^,t). As shown below, the spectral correlation function 
P{^,t) turns out to be a quantity that can be measured directly with a dedicated experiment 
from which the single emitter linewidth can be extracted. 



2.2. Measurement setup 



The experiment involves a setup previously introduced to observe the spectral dynamics of an 
isolated emitter at high temporal resolution fSl. The emission is sent to a Michelson interfer- 
ometer with an arm continuously moving back and forth (at velocity V) over a range of several 
fringes around an optical path difference 5. The intensities 4(0:4(0 the outputs of the 
interferometer oscillate as the Fourier transform of the emission spectrum 

/a,b(0 = L^,(0[1 ± I ^/((O,0cos(aj5(0/c)dco], (3) 

where 5(0 is the instantaneous optical path difference between the arms (5(0 = 5). A pair 
of avalanche photodiodes detects these intensities and a photon-counting board computes their 
cross-correlation function ^x (t) : 



WW+I), (4, 

Ia{t) 4(f + T) 



where . . . indicates time-averaging over acquisition time of the intensity correlation. 

Assuming the scanning speed V is set low enough to ensure that fringes oscillate on the 
photodiodes with a temporal periodicity c/2(0qV larger than the timescales T under investiga- 
tions, the time-averaged intensity cross-correlation function g'^ (t) measured at the output of 
the scanning interferometer decomposes as 

gX(T)=r^(T)+/"S''=(T), (5) 

with : 

r^(^) = ^[\-\vY\p^''\Q]s,^ (6) 

^^^H^) = (^gP'(T)-l+^)x (l-^FT[//-g"=(C,T)]5/,) (7) 

where /?'="='(C) =P""'(C)//^'™'(C)dC and p^'"g'<=(C,T) =P^"gi'=(C,T)//P''"gi'=(C,T)dC denote 
the normalized spectral correlation functions of the ensemble spectrum and single-emitter spec- 
trum respectively. 

This result simplifies in a few important cases. For an isolated emitter (A^ = 1), the cross- 
correlation function ^x (t) reduces to its single-emitter component, i.e. (t) ~ ^^'"^''^(t). As- 
suming the emission obeys Poissonian statistics (i.e. ^^''(t) = 1), the intensity cross-correlation 
function reads : 

1. 



(T) = 1--FT /"g'^(C,T) 



o/c 



in agreement with previous theoretical findings fS). In this case, the cross-correlation function 
g^i'c) depends on the value in 5/c of the Fourier transform (in Q of the single emitter spec- 
tral correlation function p'*'"8''^(i^, t). Hence, the dynamics of the single emitter linewidth (as 
encoded in p'*'"gi<^(C,T)) can be reconstructed through the successive accumulation of cross- 
correlation functions ^(t) at various optical path differences 5. 

Our main point is that the single-emitter component ^^'"^''^(t) in the cross-correlation func- 
tion gX (t) does not vanish even when the number of emitters involved in the experiment be- 
comes very large (A^ ^ 1), provided photoemission occurs with non-Poissonian statistics (i.e. 
^(^'(t) 7^ 1). Indeed, Equations|5]|2]then rewrite as 



^ (t) - g^^\^) - ^FT [p'^\Q + (g(2)(^) _ l)^single(^^ ^) 



b/c 
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Fig. 3. The standard FCS intensity correlation function, g^-\r). Our method works in 
situations like FCS, where the ensemble emission exhibits an intensity correlation func- 
tion ^'■"^(t) 7^ 1 at short timescales. In these cases, the single emitter spectral correlation, 
psingle^^^ t), is weighted by g'^' (t) — 1 and can be separated from the background ensem- 
ble spectral correlation, p™^{Q. 

where the normalized cross-correlation function (t) contains contributions from both the 
inhomogeneous ensemble and single-emitter spectra, with magnitudes of order 1 and (t) — 
1 respectively (Figure 3). Information on the single emitter spectrum will therefore survive 
ensemble averaging when emitters, for example, radiate intermittently, as seen on blinking or 
transiently-excited emitters - for which ^^^'(t — > 0) = 1 + 1/n, where n is the time-averaged 
number of emitters significantly contributing to the ensemble spectrum at any time. The result 
also holds for emitters showing complete photon antibunching (i.e. individually behaving as 
single-photon sources), for which g^^\T — > 0) = 1 — 1/n. 

We also note that recording cross-correlation functions (t) over long durations raises the 
signal-to-noise ratio, thus improving the measurement of the homogeneous linewidth - a result 
in complete contrast with standard spectroscopy, where longer accumulation times yield higher 
signal-to-noise ratios at the expense of greater inhomogeneous broadening. 

3. Application: unveiling spectral fluctuations in Fluorescence Correlation Spectroscopy 

Fluorescence Correlation Spectroscopy (FCS) provides a well-defined framework to illus- 
trate these findings and their significance. We choose two examples to simulate. In our first 
example, we model an ensemble comprised of single emitters with differing center frequen- 
cies, 0);, and fixed doublet spectrum consisting of two delta functions at CO,- — ^ ft), + 
As shown below, the underlying doublet is easily resolved from the broad ensemble spectrum 
with our method, despite the lack of evidence for a doublet in the ensemble emission spectrum 
(Figure 4). In our second example, we allow the center frequency of each doublet to fluctuate 
in time with a frequency a),(r) and demonstrate our ability to observe these spectral dynamics 
(Figure 5). 

3.1. Numerical methods 

We modeled FCS experiments by simulating photodetection times and emission wavelengths 
for spherical emitters of 2 nm radius freely diffusing in water at room temperature (diffusion 
coefficient D = 100 /im^/s) and excited by a tightly focused beam forming a spherical Gaussian 
excitation spot of width wq = 200 nm. The concentration of the emitters in the solution was 
adjusted so that a number of emitters « = 10 were found in the excitation volume at any time. 




Detection of the fluorescence from the excitation spot was set to a total photodetection rate of 
7=10^ counts/s. 

FCS simulations were performed by generating three-dimensional Brownian motion trajec- 
tories for emitters diffusing in a finite, cubic-shaped, open volume simulation box centered 
on the excitation spot. To do so, the diffusion trajectories of particles were first computed in 
unbounded free-space, and then subsequently put into the bounded simulation box by consider- 
ing the latter as the unit cell of a three-dimensional tiling with periodic boundary conditions. A 
concentration of n = 10 emitters under the laser spot was reached for a total number of emitters 
A'^ = 10^ in the simulation box. Once the single-emitter trajectories were computed, each of the 

emitters was assigned a center frequency o,- drawn from the underlying Gaussian ensemble 
distribution of width A centered on coq. Open volume conditions were enforced by redrawing 
the wavelength of an emitter from this distribution every time it reached a boundary of the 
simulation box. 

Computations of Brownian trajectories can be highly demanding when high temporal and 
spatial resolution are required, as is the case in FCS. However, the problem simplifies by noting 
that the position of an emitter needs only to be determined when it creates a photodetection 
event. The positions of each emitter were therefore first computed assuming uniform excitation 
over the simulation box, i.e. at times separated by intervals distributed with Poissonian statistics. 
Non-uniform excitation over the simulation box was then taken into account by filtering these 
photodetection times with a survival probability p = exp(— r^/2wQ) given by the Gaussian 
excitation profile at the emitter's location r. Emitters far from the excitation spot at a given 
time do not radiate, and so do not contribute to the FCS signal, making the computation of 
their trajectory at that time unnecessary. The width of the simulation box (4 /xm) was therefore 
kept minimal, yet large enough compared to the excitation spot size wq = 200 nm to keep finite 
simulation box effects negligible. 

Implemented in C (Anjuta 2.4.1) on a personal computer (1 GHz CPU, Linux), the above 
procedure typically required a few hours to produce long streams of photodetection events 
consisting of more than 10^ photodetection times with their associated detected wavelength for 
emitters freely diffusing in a liquid environment under focused laser excitation. The validity 
of our approach was confirmed by excellent agreement between the simulated and theoretical 
intensity correlation functions ^'^'(t) expected from FCS experiments on spherical shaped 
emitters in water (Figure 3). 

3.2. Line shape of single emitters with static spectra revealed despite ensemble broadening 

In our first example, ensemble emission was centered at a wavelength Aq = 600 nm, with in- 
dividual transition frequencies Oi , ... distributed around Aq with Gaussian statistics over 
a FWHM range A = 6 nm as expected from slight but significant inhomogeneous broadening. 
Each single emitter spectrum consisted of a doublet of monochromatic lines separated by a 
spectral width £2=1 nm. 

If measured with a conventional FCS detection setup, the correlation function §'^^(t) of 
the sample peaks at short timescales T ^ Td (Figure 3), accounting for the fact that the total 
detected intensity I{t) fluctuates as emitters continuously enter and exit the excitation volume 
with an average diffusion time To = w^/AD = 100 /is. Here, an emitter therefore diffuses out of 
the spot within a duration Td comparable to the average delay «// = 100 /xs between its detected 
photons, hence contributing to the total fluorescence signal by a few photons at most. Under 
such conditions, the spectrum observed in standard spectroscopy reduces to its inhomogeneous 
component - namely a broad Gaussian line of width A £2. The doublet in the single emitter 
spectrum is not detected (Figure 4a). 

We then add the Michelson interferometer in the photodetection path and scan the interfer- 




Fig. 4. Standard spectroscopy versus spectral correlation measurement, a) Emission spec- 
trum as measured by standard spectroscopy after an acquisition time of 1 s. The corre- 
sponding lineshape coincides with the average ensemble spectrum, {S(co,t)), and shows 
no evidence of the underlying single-emitter doublets Si{co,t). b) Using our method, the 
spectral correlation of the underlying doublet, p'*'"s'''(^, t), is easily seen on top of a broad 
ensemble pedestal, p^"'^{Q. The amplitude of t) is determined by the intensity 

correlation function g'^' (t) of the sample emission. 



ometer continuously over 10 fringes around various optical path differences 5. The scanning 
speed V is set to 5 fringes/s to access the spectral dynamics at all timescales T < c/IcoqV ~ 200 
ms. Each cross-correlation function (t) is measured by accumulating photons over 1 minute. 

The cross-correlation function ^^(t) now shows a strong dependence on the optical path 
difference 8 where it was recorded, which directly provides the normalized spectral correlation 
function p(C,t) = p™'(C) + (s'-^Ht:) - 1)p''"^'''(C,t) by taking the inverse Fourier transform 
in Eq|9] Photons separated by durations T ^ can not be spectrally correlated, since the 
population of emitters in the focal spot undergoes complete renewal over timescales T ~ To, 
accounting for the fact that t) then coincides with p™'^{Q. On timescales shorter than To, 
emitters generally do not have time to diffuse out of the excitation spot. Photons separated by 
durations T < Tq have a non-zero probability ^'^'(t) — 1 of being from the same emitter and 
therefore containing information from a single emitter It is at these timescales that the underly- 
ing doublet is revealed, as seen from the triplet of lines of relative amplitude {1/4, 1 /2, 1 /4} at 
frequencies { Q., 0, produced by the autocorrelation of a doublet of width Q. (Figure 4b). 

3.3. Line shape of single emitters with a dynamic spectrum 

The method also yields results on single emitters with a dynamic spectrum. Simulations were 
performed for emitters with a doublet spectrum of width £1=1 nm undergoing both static 
broadening over a FWHM range A = 6 nm and individual dynamic Gaussian spectral fluctua- 
tions with a correlation time Tc = 100 /xs over a broad FWHM spectral range C7 = 3 nm around 
their center frequencies CO,. All other parameters (emitter radius, scanning speed, detection in- 
tensity etc.) were left unchanged from previous section. 

In this case, the time-averaged spectrum is a broad Gaussian, while the time-resolved spec- 
trum is a doublet of separation H = 1 nm. Here again, standard spectroscopy cannot provide the 
single emitter linewidth, as spectral broadening caused by diffusion under the excitation spot 
and spectral diffusion of the single emitter happens at a rate Tjj ' + T^' = 20 ms"' faster than 
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Fig. 5. Intensity cross-correlations at thie ouputs of the scanning interferometer for a pop- 
ulation of emitters undergoing both static and dynamic spectral broadening, a) (t) as 
measured at different interferometer positions 5. The cross-correlation functions (t) dif- 
fer from the standard FCS intensity correlation function g^~\T:) at long timescales T 2> T/j 
due to the ensemble spectral correlation ij) and at short timescales T <C To due to the 
time-dependent single-emitter spectral correlation function f). b) Decay of the 

(t) with 5 at short timescale T = 100 ns (insert). The corresponding Fourier transform 
of g(2)(t:) -gX(T) for increasing values of T reveals the time dependent single emitter 
spectral correlation, p^'"s''^(i^, t). 



the average single-emitter photodetection rate I/n= 10 ms^^ 

Figure 5a shows (t) for several different interferometer positions 8. The corresponding 
patterns in the correlation functions §^ (t) can be understood as follow. For optical path differ- 
ences 8 comparable to the ensemble coherence length A = A/c = 40 jim, the output intensi- 
ties are strongly modulated by the interference pattern of the ensemble spectrum, correspond- 
ingly producing strong intensity anticorrelations between the interferometer outputs at every 
timescale T (5 = 44 /im. Figure 5a). If we now increase the optical path difference 8 until 
the ensemble coherence length A is exceeded, fringes emanating from the ensemble spectrum 
completely vanish (8 = 348 jim. Figure 5a). In this regime, distortions from the coiTelation 
function seen in standard FCS (and in Figure 3) nonetheless persist due to single-emitter in- 
terference phenomena. Indeed, when plotting ^^(t) as a function of 8, a beatnote is then ev- 
idenced at very short timescales T < 100 ns, showing that the doublet is resolved over delays 
T <C min(Tc, Td), in agreement with our previous theoretical analysis (Figure 5b, insert). Finally, 
for very large optical path differences, all interference phenomena vanish, and the standard FCS 
correlation function is actually recovered (8 = 1260 jim. Figure 5a). 

Calculating the inverse Fourier transform of g^^^ (t) — (x) for different values of T clearly 
shows the autocorrelation of a temporally evolving doublet - namely a triplet of intensities 
{1 /4, 1 /2, 1 /4} in {—£2,0, +£2} superimposed on the broad autocorrelation of the ensemble 
spectrum (Figure 5b). 

4. Conclusion 

Our approach overcomes both ensemble and temporal averaging effects in large populations of 
single emitters to provide the linewidth of a single emitter, even if many emitters are detected 
simultaneously, with each of them contributing only a few photons to the ensemble spectrum. 



We demonstrated our approach with simulations in conjunction with FCS, showing the abil- 
ity to extract a single emitter line shape from an inhomogeneously broadened ensemble. We 
then made that line shape time dependent and were able to observe the spectral dynamics with 
high temporal and spectral resolution. 

No assumption was made as to the nature of the excitation beam. Illustrated here under 
continuous excitation, our approach applies to emitters excited by a broadband lamp or a 
monochromatic laser. Pulsed excitation is also possible, particularly for the exploration of spec- 
tral correlations (e.g. multi-excitonic spectral lines) occurring on timescales shorter than the 
excited state lifetime of the emitters. 
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